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We apply the covariant derivative expansion of the Coleman-Weinberg potential to the sfermion 
sector in the minimal supersymmetric standard model, matching it to the relevant dimension-6 
operators in the standard model effective field theory at one-loop level. Emphasis is paid to nonde¬ 
generate large soft supersymmetry breaking mass squares, and the most general analytical Wilson co¬ 
efficients are obtained for all pure bosonic dimension-6 operators. In addition to the non-logarithmic 
contributions, they generally have another logarithmic contributions. Various numerical results are 
shown, in particular the constraints in the large Xt branch reproducing the 125 GeV Higgs mass 
can be pushed to high values to almost completely probe the low stop mass region at the future 
FCC-ee experiment, even given the Higgs mass calculation uncertainty. 


I. INTRODUCTION 

The first task of next generation colliders such as the ILC, the CEPC and the FCC-ee is precision measurement, of 
Higgs physics as well as Z-pole physics and tri-gauge boson physics and so on. If there is new physics, it is expected 
to show up first as corrections to the standard model (SM) processes, subject to such precision measurements. The 
corrections can be naturally sorted by dimension-6 operators in the SM effective field theory (EFT), which has a 
total of 59 independent operators for one family of fermions as a complete basis [l[. Precision measurements can be 
translated into a set of constraints on (part of) the operators. 

The operator level fitting process is (UV) model independent. On the other hand, in literature there are many 
popular new physics models with different motivations and merits. Previously the matching of new physics to the 
model independent operator constraints are mostly electroweak precision test (EWPT) and Higgs cross sections/decay 
branching ratios. For the most important subset of pure bosonic operators and at one-loop level which is usually leading 
with occasional exceptions, a method called the covariant derivative expansion (CDF) [^, Q has greatly facilitate 
the matching procedure, with advantage of being complete for all operators simultaneously, model independence in 
computation process, and giving analytical results. 

It is an industry to scan for UV models and match them to the SM EFT , but before that some generalizations 
of the CDE technique need to be made. In addition to the generalization to fermionic degree of freedom (DOF) 
which has extra gamma matrices related contributions @ , for realistic model parameters the original assumption 
of degeneracy of large scales is too simplified to be dropped. Since the CDE algebra has already been formulated 
in matrix basis, it’s not a hard generalization to use nondegenerate large scales, and use heavily the integration of 
Feynman parameterization at the textbook level after taking traces. 

Here we do sample calculations for the sfermion sector of the minimal supersymmetric standard model (MSSM). 
We expand the Coleman-Weinberg (CW) potential with covariant derivatives, matching it to a maximal set of bosonic 
operators, getting their Wilson coefficients. In addition to the nondegeneracy of soft mass squares we keep the small 
bottom/tau Yukawa couplings, and the differences between the squark sector and the slepton sector are accounted by 
different SU{3)c representations and U{1)y hypercharges. The analytical results at one-loop level should be complete 
and realistic to use for various purposes. 

A lot of works have been done in the indirect precision constraints of the squark sector of the MSSM, including the 
one-loop level Feynman diagram calculation. To provide some new useful numerics here we consider the MSSM Higgs 
mass constraint in detail, in particular we solve the mixing term Xt for the required SM like Higgs mass. There are in 
general two solutions for a required SM like Higgs mass, and we find the precision constraints are especially effective 
for the large Xt branch. 

This paper is organized as follows. In Section im we briefly review the dimension-6 operator basis and the bosonic 
CDE formulism, then we discuss the structure needed for large mass nondegeneracy. Next in Section Hill we show the 
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sfermion sector input to the CDE formulism, and then the main analytical results. Section llVI is devoted to various 
numerical applications of the analytical results. We conclude in Section [V] At last, Appendix provides a list of 
useful auxiliary formulas for the integration of Feynman parameters. 


II. FORMULISM 
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TABLE I: Independent CP-even diniension-6 operators composed of only the Higgs and gauge boson fields that are relevant to 
the analysis in this work. 

The operator basis we use is listed in Table HI which includes only Higgs and electroweak gauge boson fields. The 
basis is natural in sense that it directly captures all the possibilities generated by the CDE formulism, before removing 
any redundant operators. 

The bosonic CDE formulism given in @ can be summarized as the dimension-6 terms generated by the one-loop 
integration 
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In Eq. ([ T ]) ub is the DOF for each entry, being 1 for a real bosonic DOF and 2 for a complex bosonic DOF. 
The dimension-2 auxiliary number u is introduced as a trick to regularize orders of commutators of ^s acting on 
(p^ — M^)s 0- A Wick rotation is performed and subscript “E” indicates Euclidean. The CDE assumes in the CW 
potential the double derivatives (indicated by ') to the beyond SM fields acting on the Lagrangian potential terms 
can be decomposed as V" = -|- 5V", where is some large constant squared mass term which is irrelevant to 

Higgs vacuum expectation value (VEV), and 5V" on the other hand captures the spacetime dependent part (namely 
the SM Higgs) of potential terms, with not only the physical Higgs boson but the Nambu-Goldstone modes before 
electroweak symmetry breaking. At last ~ generally indicates a Baker-Campbell-Hausdorff expansion with covariant 
derivatives and ^s, which can be understand as introducing spacetime dependence in the momentum space @ (and 

the covariant generalization in gauge field d ^ D = d — igA is given in [^, [l^ , in which the G terms arise). 

Note that the 5V" and G are all matrices, with each entry coupling to a gauge DOF of new particle to be integrated 
out. So is the large spacetime independent mass term which in all known models is a diagonal matrix (if is 
not diagonal the new physics has intrinsic mixing irrelevant to the SM, such a model seems difficult to get motivated). 
This directly shows the way of generalizing the CDE formulism with nondegenerate large mass parameters, namely 
treating l/(p|; -I- -|- m) as a diagonal matrix (p|; -|- -|- u)“^ and doing the multiplication with full respect to their 

matrix nature and non-commutativity with the {6V" + G)s. The ^ action on any -|- + u)~^ factor (as well as 

the factor p in the first term of G) will give a commutative p number on the numerator and introduce no ambiguity. 
Eventually the overall trace operation reduces the matrix to numbers. 

The later jobs such as the textbook trick of Feynman parameterization for each term, the dimensional regulariza¬ 
tion/reduction integration over f dfipE and the integration over J du (with the the boundary u ^ oo dropped as MS 
subtraction) are quite tedious, but possible especially with a symbolic calculation tool such as the Mathematica. 
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III. THE SFERMION SECTOR 

In the following we do the CDE matching for the squark sector explicitly. Without specifying the Yq,Yt, Y^ values 
the slepton sector results can be correspondingly obtained, though one should keep in mind that the cqg vanishes 
and every other Wilson coefficient should be divided by a factor of 3 as the SU (3)c color multiplicity. 

The CW potential matrix in the basis of (ti, 6/?) has large spacetime independent masses (suppressing the 

SU{3)c components) 
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and the spacetime dependent terms can be written in a decomposition SV" = SVp + SVp + , with 
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which count the supersymmetric F-term, D-term and the trilinear X-term contributions respectively. Here the Yukawa 
couplings are defined as their SM values such as yt = V^mtjv = yJ^®®'^sin/3 with v = 246 GeV, and we have ignored 
any possible CP phases. Plugging them into Eq. © and collecting the dimension-6 operators as described in the 
appendix of we tabulate the resulting Wilson coefficient for each operator in the following long table. 

TABLE II: The Wilson coefficients of integrating out the sfermion sector. The operators are listed in Table U and the Wilson 
coefficients are defined with dimension—2. 
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Non-logarithmic Contributions 
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TABLE II: continued 



Non-logarithmic Contributions 

Logarithmic Contributions 
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Here we have a few comments: 


• In general every Wilson coefficient has a non-logarithmic contribution and logarithmic contribution, the latter 
of which is proportional to the logarithm of the soft breaking mass ratio. This is the general way the scales 
nondegeneracy comes. In that case the denominators of both contributions contain factors of certain powers of 
the nondegeneracy, with the logarithmic one higher by one in power, consistent with the behavior of the Feynman 
parameter integration listed in the appendix. It is straightforward to check that, while each contribution is 
singular in the degeneracy limit, the combination of the two are finite and reproducing an expansion around the 
degenerated results. 


• For a hierarchy of the large scales, only the non-logarithmic contribution survives, the logarithmic contribution 
will always be more suppressed by one extra large scale square. It is straightforward to read them out, for 
example the leading contribution to ct is 
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Note that in the second hierarchy even the leading contributions are suppressed by Mx. 
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• The symmetry between the up type and the down type is manifested in the top-bottom switching t o 6, which 
by definition is yt ^ yt, Mi —>• M/,, Xt —>■ Xb,Yt —>■ Y/,, and yt —t yt, M^ —>• M^, X/, —W, Y/, —>• Yt in some cases 
that the “down” type quantities also appear in a “up” type term. The SU{2)l couplings includes a induced 
flipping between the up and down sector, since they are the same for both squarks and sleptons we do not keep 
it explicitly like the hypercharge Yt(Y/,). Their effect is to flipping the sign of c/2 —/> —g^. And note the sign flip 
of the cwB, cb and chb in the top-bottom switching. 


• In addition to the above operators, there are “universal” contributions to the pure gauge boson operators. Since 
the new particle couplings to SM gauge group are always within a gauge representation so that share the same 
large scale, for example in our sfermion case the SU{2)l gauge bosons only couple to left hand sfermions so 
that only feel the Mq{Mi), and the other SU{3)c and C/(l)y gauge bosons couplings are block diagonal in the 
m2 + SV" matrix, they do not show the above nondegeneracy property. Here we adapt the results in Q for 
completeness 
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Only the operator O 3 W is directly constrained by the tri-gauge boson precision experiments, and is generally 
expected to have a relatively low sensitivity. We do not include them in the following fit. 


• In the following numerical works we will focus on the four most stringent constraints, but completely matching 
to EFT indeed has an advantage of providing much more information. For example, following [3| the hZ'j, 
hWW, hZZ can also be calculated and constrained with the Wilson coefficients. 

• The degenerated mass case result is obtained in ii- The nondegenerate mass result for cgg, cww, Cbb and 
cwB is obtained in [^, however the latter three are in a different basis. Here it is straightforward to check that 
for the Higgs diphoton coupling the combination of terms proportional to ytXf from the latter three operators 
conspires to cancel the (M| — M~)~^ and the logarithmic term, going back to a form consistent with the Higgs 
low energy theorem. The alternative Feynman diagram calculation dates back to m- 
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IV. NUMERICAL CONSTRAINTS 


The afore calculated Wilson coefficients can be used straightforwardly to transfer the model independent constraints 
to a set of constraints of the sfermion parameters. In the following we will ignore the subdominant slepton sector 
contribution, and in the squark sector we will ignore the much smaller bottom Yukawa couplings^, and for simplicity 
we further assume the right handed soft squark mass which presumably will not make significant differences 

if relaxed. While the results in TablelUdepend on the two soft breaking parameters Mg and M^, they can be translated 
to the physical masses of the two stop squarks, through mixing determined by diagonalization of stop mixing matrix 
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Note that the mixing matrix is just the first and third rows and columns of the matrix V = + 5Vp + 5Vp + 5Vx 

in Eq. (irom , while the Higgs components are replaced by their VEVs. We will ignore any loop correction to the 
above squark mass matrix. 

There are four most stringent operator level constraints in the SM EFT. Two of them are the EWPT Peskin- 
Takeuchi T and S parameters [l3|, and the other two correspond to two channels of Higgs production/decay coupled 
to gluon pairs and photon pairs, which are generated at one-loop level in the SM even in leading order. The four 
constraints in our operator basis of Table |T] are represented 0 

S' =47ru^(4cwB + cw + c_b), (11) 

XPi,-- (47r)^ o, , , , 

where Mhgg and is the SM amplitude of the (looped) hgg and hjj couplings. Note in the Higgs precision 

experiments usually measurements are not for ®'’^® subdominant 

contributions induced by other operators beyond the interference of new physics and the SM amplitude, but for 
simplicity we will ignore all of them. The current best measurements and projected future sensitivities are listed in 
Table mil 
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Current 

0.07 [13] 

0.09 yj 

4.6% yj 


ILC (ITeV) 

0.022 [1^ 

0.016 [IJ 

3.1% [1^ 

8.5% yj 

CEPC 

0.009 [1^ 

0.014 [1^ 

1.9% [1^ 

9.1% [1^ 

FCC-ee 

0.004 [1^ 

0.007 [1'3 

1.4% [18] 

3.0% 


TABLE III: The uncertainties achieved or expected at each experiments, for the four most constraining model independent 
operators or observables. 

A lot of numerical fittings have already been performed in the literature [1, [2^ [2l|. In Fig. [T] we focus on the 
2a constraints of the T and S operators as well as the hgg coupling for comparison with [J for different choices of 
Xt, with nondegenerate XIq and substituted by the physical stop masses Mi^ and The Xt introduces some 
splitting between the two stop masses, so for nonzero Xt a region with Mj:^ ~ is inaccessible and shaded gray in 
the plots (also regions with very small Mi^ and large Mi^). We assume ti has larger mixing component of left handed 
stop and t 2 of right handed stop, so that the bottom right corner of each plot shows the constraints for Mg < 
while the Mg > case is shown in the top left corner. Since the Xt = 0 slice will not satisfy the SM like Higgs mass 
constraint of 122 GeV < Mh < 128 GeV in the MSSM anyway (except for exponentially large soft mass) we will not 
show such result, and the three plots shown are for Xt = 1 TeV, 3 TeV and respectively. 

In all the three panels the most constraining operator arises between the T and hgg. On the Mg < Mi side a 
sizeable T is generated as the most stringent one, while on the other side T is suppressed and the Mg-Mi symmetric 


^ So that the which is always multiplied by yt, will not affect. In the large tan /3 region this can be not true, but for simplicity here we 
will not consider such scenario. 
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FIG. 1: The 2a T parameter (blue curves), S parameter (cyan curves) and hgg coupling (red curves) constraints, shown for 
the left Xt = I TeV, middle Xt = 3 TeV and right Xt = panels. For each group of the T, S and hgg curves the 

current constraints, ILC expected constraints, CEPC expected constraints and FCC-ee expected constraints are denoted by 
solid curves, dashed curves, dotdashed curves and dotted curves respectively, if the constraint is strong enough to be shown on 
the plot. The shaded region is theoretically inaccessible for the Xt choice with mixing indicated by Eq. (O. Here we choose 
tan P = 20 but the tan P dependence is weak. 


hgg takes over, but the T parameter is still somehow competitive if is close to , especially for large Xt and 
more aggressive experiments. In the first panel the constraints vanish at ~ Xt or ~ Xt oi 1 TeV, that’s 
the “blind spot” discussed in [2^ [2l|. Except for the blind spot, the constraint for the heavier stops will naively 
extend to large values, that’s another interesting feature of the indirect constraint method. We have checked that the 
constraints depend very weakly with tan /3, which can be understood through the fact that it always comes in the 
cos 2/3 factor of the small D-term. 

In supersymmetric theories the SM like Higgs mass is calculable, from the D-term tree level contribution and the 
large radiative corrections mainly from the top/stop sector. For stop masses of a few TeVs, the Xt term needs to 
contribute significantly to tune the Higgs mass to the measured 125 GeV, which predicts the Xt to be a considerable 
value with roughly the same scaling with stop masses, potentially contributing significant X^ and Xf terms of 
the Wilson coefficients. So the SM like Higgs mass imposes another constraint of the parameter region, and it is 
interesting to stick to that parameter slice. Note that there is still some uncertainty in calculating the SM like Higgs 
mass, especially in the low SUSY scale region, different code may give up to 3 GeV discrepancy. We use the FeynHiggs 
2.11.2 ^ to calculate the MSSM Higgs mass (for an estimation of hierarchical stops two-loop contribution, see Eq. 
(5.3) of 2^), then solve for the Xt values which are consistent with a SM like Higgs mass of 122 ~ 128 GeV. 

In Fig. [5] and [3] we show such numerical constraints for three future experiments and a range of SM like Higgs 
masses at tan/3 = 20 and decoupling large CP odd Higgs mass. Again we use the same notations and methods 


for the two cases of Mg < Mi and Mg > Mi, and the gray shaded theoretically inaccessible region. The 
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region and yellow 2a region are determined by fitting to all the EWPT and Higgs precision observables in 
(including the hWW, hZZ and Higgs to fermion pairs expectations ignored in Table iHl)) . and four individual 2cr 
constraints corresponding to Table IHII are shown as well. For a fixed SM like Higgs mass value there are in general 
two Xt solutions as shown respectively in the two sets of plots, and the most effective constraints are seen in the large 
Xt branch. While the CEPC can probe most low stop mass (< 2 TeV) region except for two corners, the FCC-ee 
can completely cover the whole region, even with a 3 GeV SM like Higgs mass uncertainty. On the large Xt branch 
one also need to worry about the possibility of other vacuum deeper than our electroweak symmetry breaking one 
with nonzero VEV of stops scalars, which may be induced by the large non-diagonal Xt- We show the bound of 
At/,/M? + M? < of 241 in Fig. [3] for comparison. 

y *1 *2 


At last we check the tan /3 dependence, which has a large effect on the SM like Higgs mass. As can be seen from 
Fig. m smaller tan /3 will allow a smaller At on the large At branch, so that weaken the constraints compared with 
the large tan/3 case. Still we can see that even if tan/3 is as low as 5, the low stop mass (< 2 TeV) region on the large 
At branch can be almost completely probed. 
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FIG. 2: The expected constraints from future ILC (the left column), CEPC (the middle column) and FCC-ee (the right column) 
experiments (see Table Hill) for fixed SM like Higgs mass value of Mh = 122 GeV (the first row), Mh = 125 GeV (the second 
row) and = 128 GeV (the third row). The Icr and 2cr allowed region are green and yellow hatched. The individual 2cr 
constraint from T, S parameters and hgg, /177 couplings are shown for each experiment again as blue, cyan, red and magenta 
curves respectively, if strong enough to be shown. There are still theoretically inaccessible region as hatched in gray. Here we 
choose tan /3 = 20, and the small Xt solution to reproduce the required SM like Higgs mass. 


V. CONCLUSION 

In this paper as an example of the CDE generalization, we perform the one-loop integration out of the sfermion 
sector in the MSSM, with full respect to all the coupling constants and nondegeneracy of the soft mass squares, 
and matching it to a basis of dimension-6 pure bosonic operators. Analytical expressions are given for each Wilson 
coefficient, with in general non-logarithmic contributions and logarithmic contributions. 

Numerically in the language of EFT, the most constraining T parameter are taken into account in a general way, 
and comparison is made among all the most stringent operators. Assuming the SM like Higgs mass relation in the 
MSSM, the probed region for each future experiments are shown. In particular in the large Xt branch the constraints 
can be pushed to very high values, and probably rule out the low scale stop sector, by precisions provided by, e.g., 
the FCC-ee experiment. 
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FIG. 3: Same as Fig. (2] but for large Xt solution to reproduce the required SM like Higgs mass. The brown curves are the 
charge and color breaking vacuum constraint of Xt/. M? + M? < and the region below and on the left are allowed. 
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Appendix A: More on Feynman Parameterization and Its Integration 

After taking trace, the denominator of each term is a multiplication of (j)\ + Mf + u). Each term eventually 
contributing to dimension-6 operators in Table U actually falls into one of the following categories 


pI 




nEm=5(p|+^"+«)"•’ 
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FIG. 4: The expected constraints from future ILC (the left column), CEPC (the middle column) and FCC-ee (the right column) 
experiments (see Table IIII[1 for fixed SM like Higgs mass value of Mh = 125 GeV, with small Xt solution (the first row) and 
large Xt solution (the second row). Other notations are the same with Fig. [2] In the large Xt branch all the theoretically 
accessible region are consistent with the charge and color breaking vacuum constraint of Vt/./M? + M~ < y/S. 

y ^ 1 ^2 


_ 

X,X,XkXi 

X^X.XkXiXmXr, 

Uj:n,=7iPl + + 


pIX,X, PeX,X, 

nEn,=6(PE + ’ nEn,=7(p|; + ’ 

plX,X,XkXi 

nEni=7(p|+^f +«)”’’ 

(Al) 


Here Xi generally indicates some dimension-1 quantity, which could be the Xt^Xf, terms in our sfermion sector, or 
the large scale itself in the vector-like fermion model. 

To see this is the case we should at first go through the following calculation steps. The standard textbook trick of 
Feynman parameterization allows the loop integration over f d'^pE, with dimensional regularization/reduction which 
makes no difference at one loop level. Then the integration over J du can always be performed, with the the boundary 
u —>■ oo dropped as MS subtraction. The first integration increase the dimension by 4 and the second integration by 
2, and all the above categories eventually give correct dimension of —2, as required for the dimension-6 operators. 
Note that each pE on the numerator corresponds to a ^ and so that a covariant derivative in the 5V" + G, checking 
with all possible combinatorics for each operator in Table H] (see the appendix of @), the only term arise from the 
Od calculation^ and no term falls into the category of p|;/„_g(p|; + Mf -|-u)”b 


^ The field strength related term such as in the combinatoric of always cancels due to the antisymmetric 

property of 
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After finishing the loop integration and f du, the Feynman parameter integration we actually need to do is 




2)712 


(E^i - 1)! 

Uirii - 1)! 


r (n_^ 


n 


„ni—l 


iE Mfx.y 


D+H ni-3 ■ 


(A2) 


Here we do not keep the factors from loop integration, and the Up and ni, n 2 , ■ • • corresponds to categories in Eq. IMl. 
In our model of sfermion there are at most three different large mass scales. For the case of two large mass squares 
we list them all 


F((p2)0, (M2)2, (M|)2) 

F((p2)0, (Mlf, (M|)2) 

F((p2)0, (m2)4, (m|)2) 
F((p2)0, (Mlf, (M|)3) 

F((p2)0, (M2)^ (M|)2) 


3(Mf-3M|) 3M| 

2(Mf - M|)2 (Mi^ - M|)3 
3(Mf + M|) 6MfM| Mf 

{Ml - M|)2 “ (Mf - M|)3 

2{Ml -5MlMl -2M^) 12M^ Ml 

M2(M2 - M2)3 + (M2 - M2)4 ^ ’ 

6(Mf + 5M|) 12Ml{2Ml + Ml) Ml 
{Ml - Ml)3 (Mf - M|)4 
5(Mf - 7MfMl - TMlMi + Ml) 30A/| 

2Ml{Ml - Ml)^ (Af2 _ ^2)5 

10(Mf + lOMi^ATl + M|) 60M|(iHf + M|) iWf 
Ml{Ml^^Miy (Afi^ _ ^2)5 

90(Mi2 + Ml) 30(Mf + AAffAfl + M|) Ml 
(iWf _ JVf2)4 + (Afi^ - Af|)5 ’ 

3Mf - 27MfMl - MMfMl + 13MfiW| - 2M| 60M| Ml 

Aff - Af|)5 (Af2 _ Af2)6 

5(3Mf + 47MfMl + llMlMl - Ml) 60M|(2Af? + 3Ml) , M? 

- M/(Mf-M|)^ -- (Mf-M|)» ‘"Ml' 

20(10M4 + IQMlMl + Ml) 60{Ml + QMlMl + 3Ml) Ml 
Ml{Ml - Ml)^ {Ml - Ml)^ ^ Ml ’ 


(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

(AlO) 

(All) 

(A12) 


F{{pY,{Ml)\{Ml)y 

F{{pY,{Ml)y{Ml)^) 

F{{py,{Ml)y{Ml)y 

F{{p^)\{Ml)\{Ml)^) 

F{{p^)\{Ml)y{Ml)y 


2{2Ml - 7MlMl + 11M|) 4Ml Ml 

3{Ml - Ml)^ {Ml - Ml)* '^Ml' 

2{Mf -5MlMl -2Ml) 12MfAf| Ml 
{Ml - Ml)^ {Ml - Ml)^ ^Ml’ 

5{Mf - TiMfMl + liMlMl + 3M|) 20M| Ml 

3Ml{Ml - Ml)^ {Ml -Ml)^ ^ 

10(Mf - SMlMl - 17Ml) 20M|(3M? + Ml) , Ml 

3(M2 - Af|)4 {Ml -Ml)5 

10(Mf + lOMfikTl + Af|) 60MlMl{Ml + Ml) Ml 

{Ml - Ml)^ {Ml - Mi)^ ^Ml’ 


(A13) 

(A14) 

(A15) 

(A16) 

(A17) 


F{{p^)y{Ml)y{Ml)y 

F{{p^)y {Ml)\ {Ml)^) 
F{{p^)\ {Ml)\ {Ml)^) 
F{{p^)\{Ml)y{MlY) 


5(3Mf - UMlMl + 23MlMl - 25M|) 5M| Ml 

12(A42 - Af|)4 {Ml -Ml)^ ^Ml’ 

5{Mf - bMfMl + ISMlMl + 3Ml) 20MlMl , Ml 

- 3(M? - M'ir - 

5{Mf - 7MfM^ - 7M}Mi + A/|) 30A/?Af| , A/? 

- 2{Ml-Ml)^ - {Ml-MD^ ^^Wr 

3Mf -17Mf Ml + Amf Ml -77Ml Ml -12MI 30M| Ml 

2Ml{Ml - Ml)^ {Ml - Ml)^ A^’ 


(A18) 

(A19) 

(A20) 

(A21) 










































F((/)^ (M2)4, (M|)3) 


5(Mf-7MfM| + 29Mi2M| + 37M|) 30M|(4Mi2 + M|) Mf 

2{Mf - M|)5 (M^ - M|)6 

5(Mf - - 47MfM| - 3M|) 60MfM|(3Mf + 2M|) 

(Mf - M|)5 (M^ - M|)6 


The three different large mass scales case can be worked out similarly. 
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